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Abstract: In this paper, we investigate the existence of Lifshitz solutions in Love- 
lock gravity, both in vacuum and in the presence of a massive vector field. We show 
that the Lovelock terms can support the Lifshitz solution provided the constants of 
the theory are suitably chosen. We obtain an exact black hole solution with Lifshitz 
asymptotics of any scaling parameter z in both Gauss-Bonnet and in pure 3rd or- 
der Lovelock gravity. If matter is added in the form of a massive vector field, we 
also show that Lifshitz solutions in Lovelock gravity exist; these can be regarded 
as corrections to Einstein gravity coupled to this form of matter. For this form of 
matter we numerically obtain a broad range of charged black hole solutions with 
Lifshitz asymptotics, for either sign of the cosmological constant. We find that these 
asymptotic Lifshitz solutions are more sensitive to corrections induced by Lovelock 
gravity than are their asymptotic AdS counterparts. We also consider the thermo- 
dynamics of the black hole solutions and show that the temperature of large black 
holes with curved horizons is proportional to Tq where z is the critical exponent; this 
relationship holds for black branes of any size. As is the case for asymptotic AdS 
black holes, we find that an extreme black hole exists only for the case of horizons 
with negative curvature. We also find that these Lovelock-Lifshitz black holes have 
no unstable phase, in contrast to the Lovelock-AdS case. We also present a class of 
rotating Lovelock-Lifshitz black holes with Ricci-flat horizons. 
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1. Introduction 

Candidates for gravity duals of non-relativistic scale invariant theories have recently 
attracted a great deal of attention for several reasons. One is that some condensed 
matter systems realized in laboratories at their critical points are described by non- 
relativistic conformal field theories. Furthermore, the AdS/CFT correspondence ^ 
[^ describing the duality between strongly coupled conformal field theory and gravity 
continues to find more applications in other branches of physics such as QCD quark- 
gluon plasmas 0], atomic physics, and condensed matter physics ||, ||, ^ ||, ^. 
This has led to a study of gravity-gauge duality in a much broader context than 
its original AdS/CFT formulation, extending to non-relativistic and Lifshitz field 
theories, with the aim of gaining more knowledge of the strong coupling behavior 
of these other physical systems. Non-relativistic conformal symmetry contains the 
scaling symmetry 

t XH, x ^ Ax (1.1) 
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where z is the dynamical exponent. This transformation exhibits anisotropic scale 
invariant behavior. Actually, for condensed matter applications, one typically is in- 
terested in considering anisotropy between different spatial dimensions. For z = 1 
this scaling symmetry is the familiar conformal symmetry. Such a non-relativistic 
scale invariance ( |1.1| ) can be exhibited by either a Galilean-invariant theory or a 
Lifshitz-invariant theory. From a holographic point of view, this suggests the follow- 
ing (asymptotic) form for the spacetime metric 



2zj.2 I '^^^ I 2 T 2 

-r dt H TT- + r ctx 



that obeys the scale invariance 

t A^'t, r -)■ A~V, X Ax 



;i.2) 



:i.3) 



noted previously in a braneworld context |]T0 |. 

A four-dimensional anisotropic scale invariant background using an action in- 
volving a two form and a three form field with a Chern- Simons coupling 



d X\ 
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can be engineered to yield solutions with this asymptotic behavior |TT]], where H = 
dB and F = dA. Such an action was argued to be rather generic in string theory, 
although no explicit brane realization or embedding into ten-dimensions was given. 
For these matter fields, lots of effort has been expended in extending this solution 
to the case of asymptotic Lifshitz solutions. One of the first analytic examples 
was reported in Ref. [0 for a sort of higher- dimensional dilaton gravity without 
restricting the value of the dynamical exponent z. An exact topological black hole 
solution with hyperbolic horizon which happens to be asymptotically Lifshitz with 
z = 2 was found in |T3|; further solutions with z = 4 and with spherical topology 
were subsequently obtained However in general such asymptotic Lifshitz black 
holes must be investigated numerically |T3|, |T6[. Other possibilities for the 

matter needed to support such a background have been investigated [|T^. Asymptotic 
Lifshitz solutions in the vacuum of higher-derivative gravity theories (with curvature- 
squared terms in the action) have been investigated 0, [l^, |19|, the higher 
curvature terms with suitable coupling constant play the role of the desired matter. 

Recently gravity theories including higher powers of the curvature, particularly 
curvature-cubed interactions, have attracted increased attention 112111. This is be- 



cause, in the context of AdS/CFT correspondence, corrections from higher powers 
of the curvature must be considered on the gravity side of the correspondence in 
order to scrutinize CFTs with different values for their central charges. Amongst 
the cornucopia of higher-curvature gravity theories. Lovelock gravity theories play a 
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special role in that the number of metric derivatives in any field equation is never 
larger than 2. Third-order Lovelock gravity is supersymmetric, and therefore one 



can define superconformal field theories via the AdS/CFT correspondence ||2^, ^ 
Furthermore quasi-topological gravity including curvature-cubed interactions, while 
not supersymmetric, can be considered to be dual to some non-supersymmetric but 
conformal gauge theory in the limit of a large number of colours [El 



In this paper, we consider the existence of Lifshitz solutions in third order Love- 
lock gravity both in vacuum and in the presence of a massive vector field. Since 
the higher curvature terms appear to play the role of some kind of matter field, it is 
natural to ask whether they can support a Lifshitz solution in vacuum or not. We 
find that the answer is yes, albeit under restricted circumstances. We also search 
for asymptotic Lifshitz black holes in the presence of a massive vector field, whose 
action is given via a dualization of the action (|1.4| ). The solutions we find - both 
analytically and numerically - can be regarded as higher-curvature modifications to 



those obtained from Einsteinian gravity coupled to matter [13, 14, 16 



The outline of our paper is as follows. We give a brief review of the field equations 
of third order Lovelock gravity in the presence of massive vector field in Sec. ^. In 
Sec. ^ we present the {n + l)-dimensional exact Lifshitz solutions in vacuum and 
in the presence of a massive vector field. In Sec. ^ we obtain the series solutions 
to the field equations near the horizon, while the series solutions at large r will be 
given the Appendix. We then obtain numerical solutions to these equations. The 
thermodynamics of these Lovelock-Lifshitz black holes will be given in Sec. ^ Section 
1^ will be devoted to the rotating Lovelock-Lifshitz solutions. We finish our paper 
with some concluding remarks. 



2. Field equations 

The fundamental assumptions in standard general relativity are the requirements of 
general covariance and 2nd-order differential field equations for the metric. Based on 
the same principles, the Lovelock Lagrangian is a very general Lagrangian in classical 
gravity that produces second order field equations for the metric for arbitrary powers 
of the curvature [^. The action of third order Lovelock gravity in the presence of 



an Abelian massive vector field may be written as 

1 = J r+'xV^ (^-2A + £i + a2C2 + a^Cs - ^F^uF'^'' - ]^m'A^A^^ , (2.1) 

where F^^ = d[f^A^], A is the cosmological constant, 0^2 and as are Gauss-Bonnet 
and third order Lovelock coefficients, Ci = R is the Einstein-Hilbert Lagrangian, 
^2 = Rf.u^sR'"'^^ - 4:Rf,^Rf"' + R^ is the Gauss-Bonnet La grangian, and 
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d3 I r)Tyiii/aKTy JDPT i Q DUV Dctk JDPt i cyAjyiiuaKT) JDp 

3 — Jt + In'^ tianprn, + crp^ vt^ fiK + ^4Jt'^ riaKupn-^ 

+mR^'''''^R,^^, + 2AR^''''^R„^R,, + 1QR>^''R,^R\ - 12RR^^R^, (2.2) 

is the third order Lovelock Lagrangian. We assume that the Gauss-Bonnet coefficient, 
which has the dimension of (length)^, is positive as in the heterotic string theory [^]. 

In Lovelock gravity only terms with order less than [{n + l)/2] (where [x] is the 
integer part of x) contribute to the field equations, the rest being total derivatives in 
the action. For 3rd-order Lovelock gravity we therefore consider {n + 1) -dimensional 
spacetimes with n > 6 (though in situations where we set 03 = our solutions will 
be valid for n > 4). Varying the action with respect to the metric tensor g^j_u and 
gauge field the equations of gravitation and gauge fields are 

G^il + «2Gg] + a,Gf, + kg,, = (2.3) 
d, {^F^'^) = m'^A-", (2.4) 



where 



is the energy-momentum tensor of gauge field, G^p} is just the Einstein tensor, and 
G')S, and G^pv are given as: 

G^pv 2( R^ui^T-R ^ IR^piifjRP 2,R^fjR ^ -\- RR^y) '^^29pui 

Gpu - -'i[AR RanXpR yrp ' ^R X^R rpR vpK + ^Ry R^nXpR rp 
J^CTKTpJ^up "T on ycrp^ rp^ K ' "-'^ VTK^ ap^ p 

\A p TCTK p pp AT) T<TK p pp _i_ A TJTpaK p p I O D D '^TP p 

+ 8R''uppR''aR''T - ^R'^urpR'^aRp " ^R'^^^^R" ^R^p - ARR'^^^^R^^ 

-\-A:K^''RprRup, — SR^yRrpR^p^ + ARR^pR^ p — R^Ryp) — -C^gpu- 
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The metric of an {n + l)-dimensional asymptotically Lifshitz static and radially 
symmetric spacetime may be written as: 

ds' = -—f{r)de + ^^+r'dn' (2.5) 
l^zz r^g[r) 



where the functions /(r) and g{r) must go to 1 as r goes to infinity. In Eq. ( |2.5| ) 
dQ"^ is the metric of an (n — 1) -dimensional hypersurface with constant curvature 
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{n — l){n — 2)k and volume Vn-i- We can write 



n— 1 j— 1 

ddf + E n sin^ djdOf k = 1 

i=2 j=l 

^^^i + E n sin^ Ojdef A; = -1 

i=3 j=2 J 

n-1 

E d9f k = 



(2.6) 



though it should be noted that our solutions are valid whenever dfl"^ describes any 
Einstein space. 

Using the ansatz 

At = q'^Hr) (2.7) 

for the gauge field and defining a2 = {n — 2){n — 8)0:2 and a^ = {n — 2)...{n — 5)0:3 
for convenience, the field equations (|2.3|) and ( p. 41) reduce to the system 



2r^h" - r [(In /)' - (In g)'] {rh' + z) - 2{n + z)rh' + 2{n - l)z 



2mH 



2/2 



(2.8) 



{n{n - l^g + {n - l)rV + 2A/V - {n - l){n - 2)kf} 
+ {n - l)a2l^r^{kP - r^g) {nr^g + 2rV - (n - 4)kl^} 

+ {n - l)a3{kf - r^gf {nr^g + 3rV - {n - 6)kl^} = 2/VT/, (2.9) 
rH^ {{n - l){n -2 + 2z)r^g + {n - l)r^^(ln/)' + 2A/V - (n - l){n - 2)kf} 
+ (n- 1)02/ V(A;/2 -r^c/) {{n ~ 4 + 4:z)r^g + 2r^g (In f)' - (n - 4:)kl^} 
+ {n - l)a3{kf - r^gf {{n - 6 + Gz^g + Sr^gilnf)' - (n - 6)kl^} 



2l^r^T;, 



(2.10) 



where prime denotes the derivative with respect to r and T/ and are 



4Pf 



4Pf 



[girh' + zhf + mH''h?] 
{g{rh' + zhf - mH^h^] 



3. Lifshitz Solutions 



3.1 Vacuum Solutions 

We first investigate the possibility of having (n + l)-dimensional Lifshitz solutions 

„22 ;2^^2 



i=l 



(3.1) 
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in Lovelock gravity in the absence of matter. In order to have an asymptotically 
Lifshitz solution in 3rd-order Lovelock gravity, the following constraints on the cos- 
mological constant and third order Lovelock coefficient 

^ ^ ^2/2 _ , d3 = -^(/'-2d2). (3.2) 

hold for an arbitrary value of as is easily obtained via straightforward calculation. 
Note that for = the constraints (B.2D become 



n(n — , , 

A = - and «2=2- (3-3) 

and so asymptotically Lifshitz solutions exist in the vacuum of Gauss-Bonnet gravity 
{as = 0) as well. The cosmological constant here is half that of an AdS spacetime. If 
we set A = then a Lifshitz solution exists in third order Lovelock gravity provided 

ds = 2f, da = l^. (3.4) 

Thus, one may have Lifshitz solutions in Lovelock gravity without matter, demon- 
strating that its higher curvature terms can, for the proper choice of Lovelock coef- 
ficients, have the desired effect that matter fields induce. 

The preceding analysis was for k = 0. For A; = ±1 we find for any value of z that 

g{r) = l + — (3.5) 

furnishes an exact Lifshitz solution to Lovelock gravity without matter, provided one 
of the above conditions ( p.2|) , ( p.3|) is satisfied, with the function /(r) undetermined 
by the field equations (though boundary conditions constrain it to asymptote to 
1). We can choose /(r) = g{r). This is a naked singularity for A; = 1 but is an 
asymptotically Lifshitz black hole for k = —1. The arbitrariness of /(r) is due to 
a degeneracy of the field equations; if either of ( |3.2| ) or ( |3.3| ) hold then the field 
equations are each proportional to the factor g{r) — 1 — kP/r"^. 

This degeneracy of the field equations has been noted previously in 5-dimensional 
Einstein-Gauss-Bonnet gravity with a cosmological constant where (converting 



notation appropriately) condition (|3.3|) was obtained. Here we see that a more general 
degeneracy occurs in 3rd order Lovelock gravity. We shall see that this degeneracy 
is lifted when matter is present. We find that no other exact solutions to the field 
equations exist for these symmetries and asymptotic behavior. 

3.2 Matter Solutions 

Consider next the case of Lifshitz solutions in the presence of a massive gauge field 
A'^. The metric ( ^.ip and the gauge field ( |2.7D with h{r) = 1 solve the field equation 
provided 

n (n — l)z 
m' = ^—jr-. (3.6) 
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In order to have asymptotically Lifshitz solutions in 3rd-order Lovelock gravity 
in the presence of matter, the following constraint 



2 2(z-l)L4 



[{z - If + n{z - 2) + n^jL^ + n{n - l)(d2/^ - 2a^ 



^ ^ - (3.7) 
must hold for the cosmological constant and charge where we define 

L^ = l^- 2Pa2 + 3d3 (3.8) 

for simplicity and we use this definition throughout the paper. Since > we 
obtain the constraint 

7r<5 + 3^ (3.9) 

which in turn yields 

{z - l){n + z- + Sdg) + n{n - + da) ^ ^ ^ n{n - 1) A _ "3 



2/6 - 4/2 V I* 

(3.10) 

provided d2 > 0. Equation (p.lOj) shows that for ds < /^ the cosmological constant 
is negative. For ds > /^, the cosmological constant can be positive provided 

{z-l){n + z- 1)(/^ + Sda) + n{n - 1)(/^ + dg) ^ 1 

{2{z - l){z + n - 1) + n{n - ^ P ^ 2^ 2/^' 



where the last inequality comes from the condition 

Note that for Einstein gravity (d2 = da = 0) L = / and the above conditions 
reduce to those which are given in for n = 3. In Gauss-Bonnet gravity (as = 0), 
these conditions become 



2 2(z-l)(/2-2d2) 



zP 

[{z - I f + n{z - 2) + n^]{P - 2d2) + n{n - l)d2 ^ . 

A- • (3.iij 

where now P must be larger than 2d2. If d2 = /^/2, then the charge q becomes zero 
and the conditions ( |3.11| ) reduce to the conditions ( |3.3| ) as expected. More generally, 
if L = 0, then q vanishes and the conditions ( p. 71) reduce to the conditions ( p.2|) as 
expected. 

In the absence of a cosmological constant, the Lifshitz solution exists provided 



2_ 2n{n-l){z-l){a2-2P) 
~ zP[3{z - 1)2 + 3nz + n{n - 4)] ' 



„ _ (2d2 - n [{z - 1)2 + nz]+n [{n - 3)d2 - [n - 2)P] 

[3{z - If + 3nz + n{n - A)] ^ ' 

In this case the conditions ( |3.12| ) reduce to the condition (|3.4| ) if d2 = 2/2, leading 
to g = and the absence of matter. 
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4. Asymptotic Lifshitz Black holes 



In this section we seek black hole solutions in 3rd-order Lovelock gravity that are 
asymptotically Lifshitz. We therefore consider the field equations (|2.8| )- (|2?T0| ) with 
the conditions (|3.7|), but with the more general ansatz (p.5|) and (|2T 



4.1 Series solutions near the horizon 

We now consider the near-horizon behavior of such solutions. Requiring that /(r) 
and g{r) go to zero linearly, that is 

/(r) = /i {(r - ro) + ^(r - r^)' + h{r - r,f + U{r - r^f + ...}, 
g{r) = gi{r - ro) + g2{r - r^f + g^{r - tq)^ + g^ir - rof + 
H^) = /i^^ {^0 + hi{r - ro) + h2{r - ro)^ + hs{r - rof + ^{r - tq)^ + ...}, 
and inserting these expansions into the equations of motion arising from Eqs. 



]9D with the conditions ( |3.7| ), and solving for the various coefficients, we find that 
Hq = 0. This is consistent with the fact that the fiux dA should go to a constant at 
the horizon. Also, one may note that by scaling time we can adjust the constant /i 

1 /2 

by an overall multiplicative factor (note the use of in the expansion of h{r) as 
well, which is due to dt in the one-form A). We find the following constraint on the 
1st order constants 



gi = 4{3«3 [hi\z - iW + k'z{n - 1)1^] 

+2ro'/2d2 [-hi^z - l)ro^ + Pkz{n - 1)] + r^H^ [h^ {z-l)ro + z{n- 1)] 

X I { [3{z - 1)2 + 3nz + +n{n - 4)] Tq*' + k{n -l){n- 6)1^} 

-ro^fa2 { [2{z - if + 2nz + n{n - 3)] tq^ + k'^l\n - l){n - 4)} 

+roH^ [{z - If + nz + n{n - 1)] V + (n - l){n - 2)kl^\, (4.1) 



which means that not all boundary conditions are allowed. Note that here the 
coefficient hi is arbitrary, and should be chosen suitable in order to explore the 
numerical solutions. The higher order coefficients may be found easily, but their 
expressions are very lengthy and so we don't write them here for reasons of economy. 

One may also investigate the behavior of the metric functions at large r. We 
relegate this to the Appendix, where we show that the powers of 1/r for k = are 
in general non-integer. 

4.2 Numeric Solutions 

In order to find the numeric black hole solutions of the system given by Eqs. ( p.8[) - 
( p.lO|) with the conditions ( p.7|) , we define 
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and find the first derivatives of j(r), /(r) and g{r) as 

dj {n — l)zh — [z{n + z — 2)h + {n + 2z — l)rj]g {z — l)L^{zh + rj)r'^h'^ 

dr r'^g r'^fgH ' 

(4.3) 



df 



{n -l)[n + 6{z - l)]za^r^fg^ 



dr {n — \)zr^gH 

+ {n - [?,ka^{n + 4^ - 6) + (n + 4^ - 4)a2r^] /c/^ 

-{n - 1)^/^2 [3(n + 2z- 6)k'^a3 + 2{n + 2z - 4)^2^^ + {n + 2z - 2)r^] 
+zd3 {[3(z - 1)^ + 3nz + n{n - 4)]r^ + - l)(n - ^)kf'] f 
-z&ifr^ {[2{z - 1)2 + 2nz + n{n - 3)]r^ - (n - l)(n - 4)A;Y} / 
+z/V { [(2 - 1)2 + nz + n{n - \)Y + (n - l)(n - 2)kl^] f 

-{z - l)L^r^[{zh + rjfg + {n - l)zh'^]\, (4.4) 



dg 1 



|_^(^_l)^a3^6^/ + (^_l)^/V[3M3(n-2)+nd2r2]// 

A)k^a3 + 2{n - 2)^2^^ + nr^] fg 
+za3{[3{z - 1)2 + 3n2 + n(n - 4)]r^ + (n - l)(n - 6)fc/^} / 

2 {[2(z - 1)2 + 2nz + n{n - 3)]r^ -{n-l){n- 4)kY} f 
+;z/2^2 _ ly ^nz + n{n - l)]r2 + {n - l){n - 2)kl^} f 

-{z - 1)LV[(z/i + rjYg + {n - l)z/;.2]|, (4.5) 
where H is 

H = 3a3{kP - r^gf + a2pr\kP - r^g) + /V. 

Now having a system of 4 first order ordinary differential equations, one may explore 
the numeric solutions by choosing suitable initial conditions for /(ro), 5'(ro), h{rQ) 
and j{ro), where ro is the radius of horizon. This can be done by using the method 



explained in [0. First, we apply the numerical method for 2; = 1 which can be solved 
exactly, and compare the numerical solution to the exact one. Indeed, for z = 1, the 
charge vanishes, and the problem reduces to the case of black holes in third order 
Lovelock gravity with specific choices of Lovelock coefficients. In this case, we obtain 



the exact solutions of Lovelock gravity |28], which have /(r) = g{r). The diagrams 



of the functions /(r) and g{r) versus r for z = 1 with k = and k = ±1 have been 
shown in Fig ^]l|. In this figure the exact solutions have also been shown. As one 
can see from the plots, the numerically obtained functions /(r) and g{r) are equal 
and both of them are exactly lie on the exact solution within limits of numerical 
precision. We regard this as a good test of our methods. Note that for A; = 0, — 1 
and z = 1 (the solution without matter) , the metric functions increase from zero at 
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f 




0.2- 



oH 1 1 1 1 1 1 1 1 

12 3 4 

r 

Figure 4.1: /(r), g{r) and exact solution versus r for n = 6, I = 1, z = 1, 02 = -25, 
03 = .3 and rg = 0.85 for k = 1 (top), k = (middle) and k = —1 (bottom). Note that 
/(r), g{r) and the exact solution lie on each other for each k. 



r = ro to 1 at r = 00. However for k = 1 and -2 = 1, the metric functions grow to be 
larger than unity at intermediate values of r and asymptote to 1 at infinity. Now, 
we consider the solutions for 2; 7^ 1, for which matter exists, where we must solve the 
system of equations numerically. Figure |4.2| shows the function /(r) as a function of 
r for a small black hole of Einstein, Gauss-Bonnet and Lovelock gravity with z = 1 
and z = 2. As one can see from this figure, while /(r) for Einstein, Gauss-Bonnet 
and Lovelock with z = 1 (solution without matter) are almost the same, they are 
different for 2; = 2 for the same values of 02 and 03. Furthermore, for z = 2, the 
metric function /(r) rapidly increases to values larger than 1, eventually asymptoting 
to 1 as r —7- 00, while for z = 1 this function monotonically increases from zero at 
r = ro to 1 at r = 00. Thus, it appears that asymptotic Lifshitz solutions are 
more sensitive than are their z = 1 (AdS) counterparts to the corrections induced 
by Lovelock gravity. This feature also occurs for large black holes with k = —1 (see 
Fig. |4.3| ). Note that for k = —1, the function /(r) in Gauss-Bonnet gravity is larger 
than that in third order Lovelock gravity. This is due to the fact that the slope of 
the function /(r) decreases as da increases (see the terms in df /dr with the factor da 
in Eq. ( [4.5| ) with k = —1, which are negative). Figures [4.4| - [4.5| show that this fact 
also occurs for the small and large black holes with k = 0. For the case of = 1, one 
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Figure 4.2: /(r) versus r for k = —1, n = 6, I = 1, 02 = -25, 03 = .3, ro = 0.92 in 
3rd order Lovelock, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and 
blue (dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves, which are nearly identical and so lie on top of each other, correspond to z = 1. 



can see from Figs [4.6| and [4.7| that both the asymptotic AdS and Lifshitz black holes 
are different for various order of Lovelock gravities, although the difference between 
the metric function /(r) for various order of Lovelock gravity becomes more relevant 
as z increases. 

We find that the dependence on z of g{r) is not as significant as for /(r). For 
this and reasons of economy, we plot only the function g{r) in one case (see Fig. |4.8| ). 
We also find that g{r) is monotonically increasing for topological black holes with 
k = 0, whereas (as noted previously) /(r) has a maximum larger than 1 at r several 
times vq. 

We pause to comment on the relative importance of the various terms in the 
Lovelock action. While there is no essential problem in considering all terms to be 
similar in magnitude, we obtain some insight into Lovelock gravity by treating it as 
an approximation (in powers of the curvature) to a full quantum theory of gravity. 
For simplicity we shall consider solutions obtained in this paper for k = 0, which 
asymptote to the Lifshitz background ( |3.1| ). In this case the pth order Lagrangian of 
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Figure 4.3: /(r) versus r for k = —1, n = 6, / = 1, 02 = -25, 03 = .3, tq = 8.5 in 
3rd order Lovelock, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and 
blue (dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves correspond to z = 1. 

Lovelock gravity for the metric (|2.5| ) may be written as: 



where = {ap/ap)y/^Cp. The above expression for the Lagrangian is correct for 
any value of p, though in this paper we have considered only p < 3. The Lagrangian 
( [4.6|) can be evaluated on the solutions at large r by the use of the expansions of the 
metric functions which are given in the Appendix. It is a matter of calculation to 



= (-1) 



{2p - l)/2+2p-V(2p-l)(2--l) 



(n-1) 



(f)'"''|:{-""^^^^-^Vn> (4-6) 



show that the ratio CJ^j^^j CJ^ at large r is 




•p+i 



(2p- l)[n + 2(p+l)(^-l)] 
(2p+ l)[n + 2p(2- 1)]Z2 



(4.7) 



and so /^p+i < Cp/t^ . 
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Figure 4.4: /(r) versus r for k = 0, n = 6, I = 1, a2 = .25, as = .3, rg = 0.85 in 3rd 
order Lovelock, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and blue 
(dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves correspond to z = 1. 



5. Thermodynamics of black holes 



The entropy of a black hole in Lovelock gravity is |29 



S = -^Y.^akj d^~'x^Ck-i, (5.1) 



fe=i 

where the integration is done on the (n — 1) -dimensional spacelike hypersurface of 
the Killing horizon with induced metric g^^, (whose determinant is g), and Ck is the 
kth order Lovelock Lagrangian of g^u- It is a matter of calculation to show that the 
entropy of a black hole per unit volume of the horizon in third order Lovelock gravity 
is 

.^iv(i+ ^f'-;>r + ^f'-;>r ). (5.2) 

4 V - 3)rg {n - 5)r^ / ^ 

This reduces to the area law of entropy for a2 = as = 0. 
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Figure 4.5: /(r) versus r for /c = 0, n = 6, Z = 1, 02 = -25, 03 = .3, ro = 8.5 in 3rd 
order Loveloclc, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and blue 
(dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves correspond to z = 1. 

One can obtain the temperature of the event horizon by using standard Wick- 
rotation methods, yielding the result 



For the case of asymptotic Lifshitz black hole with (identified) hyperbolic hori- 
zon, the temperature of the exact solution (|3.5|) with horizon radius ro = / in vacuum 



The temperature of more general Lovelock-Lifshitz black holes can be calculated 
numerically. First, we review the thermodynamics of asymptotic AdS black holes 
(z = 1) to compare it with the asymptotic Lifshitz black holes. As one can see from 
figure the logarithm of the temperature of black holes versus the logarithm of 
the entropy for large black holes is linear, while for small black holes this occurs 
only in the case of = 0. For k = —1, one encounters with an extreme black hole 
in Einstein (blue dashed line) and Lovelock gravity (blue solid line). However, the 




(5.3) 
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Figure 4.6: /(r) versus r for k = 1, n = 6, I = 1, a2 = .25, as = .3, rg = 0.85 in 3rd 
order Lovelock, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and blue 
(dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves correspond to z = 1. 



horizon radius of the extreme black hole of Lovelock gravity is smaller than that 
of Einstein gravity. Note that the slopes of the graphs of log T versus log S for 
small black holes of Einstein and Lovelock gravity with k = 1 are negative, and 
therefore they are unstable. The unstable phase in Lovelock gravity occurs for a 
smaller radius black hole than in Einstein theory. Plotting log T versus log S for the 
case of asymptotic Lifshitz black holes with z = 2, we see from Fig. p.2| that the 
temperature of black holes in Lovelock gravity with given entropy S is smaller than 
the temperature of black holes of Einstein gravity with the same entropy. Conversely, 
at a given temperature the entropy of the Lovelock black holes - extremal and non- 
extremal - is larger than for the Einstein ones for both z = 1 and z = 2. The horizon 
radii of extreme black holes in both Einstein and Lovelock gravity for A; = —1 with 
z = 2 are smaller than their z = 1 counterparts, and for a given z smaller in Lovelock 
gravity than in Einstein gravity. 

Numerical calculations show that there is no unstable phase for z = 2 black holes 
for k = 1. We also find that for k = and the large black holes for = ±1, the 
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Figure 4.7: /(r) versus r for A; = 1, n = 6, Z = 1, 02 = -25, 03 = .3, ro = 8.5 in 3rd 
order Loveloclc, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and blue 
(dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves correspond to z = 1. 



temperature is proportional to Tq, while the entropy is proportional to Tq ^, yielding 

T oc (5.5) 

for the dependence of T on S. In order to say more about the temperature, one should 
compute the energy density of these black holes using (for example) the counterterm 
method. While some work has been done on this for asymptotic Lifshitz solutions 
of Einstein gravity [^, ^|, the corresponding formalism needs to be developed for 



Lovelock gravity, an endeavour that we hope to address in the future. 



6. Rotating Lovelock-Lifshitz solutions 

In this section we endow the Lifshitz spacetime with a global rotation. We first 
consider the Lifshitz solution with one rotation parameter. We write the metric as 

i=2 
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Figure 4.8: g{r) versus r for k = —1, n = 6, I = 1, a2 = -25, = .3, tq = 0.85 in 
3rd order Lovelock, Gauss-Bonnet and Einstein gravities red (solid), green (dotted) and 
blue (dashed), respectively. The three upper curves are /(r) for z = 2 and the three lower 
curves, which lie on each other, correspond to z = 1. 

In order to add angular momentum to the spacetime, we perform the following 
rotation boost in the t-9i plane: 

t^Et- adi 01 ^ EOi - ^t, (6.2) 

where a is the rotation parameter and S = 1 + a^/P. Substituting eq. ( |6.2|) into eq. 
( |6.1| ) we obtain 

i=2 

The transformation (|6.2|) generates a new metric if 9i is periodically identified 
since the transformation (|6.2| ) can be done locally but not globally ||32|. The periodic 



nature of 9i allows the metrics ( |6.1|) and ( p.3|) to be locally mapped into each other 
but not globally, and so they are distinct. It is a matter of straightforward calculation 
to show that the metric ( |6.3| ) is a solution of Lovelock gravity, provided one of the 
conditions (B.2|)-(B.3|) holds. Furthermore the metric (16. 31) is also a solution to the 
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Figure 5.1: logT versus logS" for n = 6, / = 1, z = 1 in Einstein and Lovelock (a2 = -25, 
as = .3) gravities with k = 1 (red), A; = (green) and k = —1 (blue), where the dashed 
and solid lines are curves in Einstein and Lovelock gravities, respectively. 



field equations in the presence of the vector field 

A = q'^ {Edt - adOi) , (6.4) 

provided the condition ( p.6|) and one of the conditions ( |3.11|) - (|3n2|) are satisfied. 

It is straightforward to generalize the metric ( p. 3D to a metric with more rotation 
parameters. The rotation group in [n + 1) dimensions is SO{n), and therefore the 
number of independent rotation parameters is [n/2]. The generalized solution with 
m < [n/2] rotation parameters can be written as 

m \ ^ 2 ™ 

ds'^ ~ ~77 ( '^^^ ~ ''^^ciidOi 1 + "^^^ [ttidt — EfdOi^ 

i=l / i=l 

2 m n—l 

-^—Y^{ad9j~ajddif + r^ ^ dO^, (6.5) 
' . . . ^ 




where S = yjl + af/P and the angular coordinates are in the range < < 27r. 

This generalization can be extended to the Lovelock-Lifshitz /c = black hole 
since the transformation (|6.2|) does not change the r-dependence of the metric func- 
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Figure 5.2: logT versus logS" for n = 6, I = 1, z = 2 in Einstein and Lovelock (a2 = -25, 
as = .3) gravities with k = 1 (red), A; = (green) and k = —1 (blue), where the dashed 
and solid lines are curves in Einstein and Lovelock gravities, respectively. 



tions and the vector field. Hence 



k \ ^ 2 



r'^q(r) P ^ 

^ ^ ' i<j i=k+l 



with the vector field 

A = q'^Kr) i^dt - Yj'idd, ) , (6.7) 



is a rotating black hole solution to the field equations provided /(r), g{r) and h{r) 
are chosen to be the functions calculated in section |4.2| for = 0. 

7. Concluding Remarks 

It is known that Einstein gravity in the presence of a massive vector field can support 
Lifshitz solutions [^. We have demonstrated here that such solutions exist in pure 
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Lovelock gravity provided the coupling parameters are chosen appropriately. The 
higher curvature terms appear to play the role of some kind of matter field, whose 
nature depends on the constants of the theory. We also found for any value of z 
an exact vacuum asymptotically Lifshitz solution ( p.5|) for = ±1, which is a black 
hole for /c = — 1, and a naked singularity for k = \. This solution exists in both 
Gauss-Bonnet and 3rd order Lovelock gravity, depending on the choice of coupling. 

We also obtained a broad class of solutions for Lovelock gravity coupled to a mas- 
sive vector field. After demonstrating that Lovelock gravity can support a Lifshitz 
solution in the presence of a massive vector field, we searched for asymptotic Lifshitz 
black holes in the presence of a massive vector field. Our numerically obtained solu- 
tions generalize those obtained in Einsteinian gravity coupled to a massive vector field 



1^, 0, 0. We found that asymptotic Lifshitz solutions (2; > 1) are more sensitive 
to the corrections induced by Lovelock gravity than are their z = \ counterparts. 

We also considered the thermodynamics of the black hole solutions. We found 
that, as in the case of asymptotically AdS black holes of Lovelock gravity, one can 
have an extreme black hole only for k = —\. However the horizon radius of both 
the extreme black holes of Einstein and Lovelock gravity for k = —1 with z = 2 
are smaller than their z = 1 counterparts. That is, the radius of the extreme black 
holes decreases as z increases both in Einstein and Lovelock gravity. Also, Lovelock 
terms decrease the radius of extremal black holes compared to their Einsteinian 
counterparts. The temperature of a Lovelock-Lifshitz black hole with entropy S is 
smaller than the temperature of a Lifshitz black hole in Einstein gravity with the same 
entropy. We also found, numerically, that as z increases the temperature of a black 
hole with a fixed horizon radius increases. Indeed, the temperature is proportional 
to Tq for black holes with zero curvature horizon and also for large black holes with 
nonzero horizon curvature. In these two cases the temperature is proportional to 

gz/{n-l)^ 

We investigated only the cases with z > 2 and found that these solutions are 
thermodynamically stable. This fact can be seen by diagrams of temperature versus 
the entropy. This feature of asymptotically Lifshitz black holes is different from the 
asymptotically AdS solutions that can have an unstable phase. 

Our results demonstrate that higher-order curvature corrections to "Lifshitz 
holography" are under control and have a sensible physical interpretation. How- 
ever there is much remaining to explore for these kinds of black holes, including the 
corrections they induce in the dual boundary theory and a more detailed study of 
the behavior for larger n and larger z. The rotating k = solutions should have 
counterparts for the k = ±1 cases. Work on these areas is in progress. 

8. Appendix 

Here, we investigate the behavior of the metric functions at large r. For the case of 
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A; = the powers of 1/r may be non- integer, and therefore we consider it separately. 
For this case, we use straightforward perturbation theory, writing 



and then finding the field equations up to the first order in e. We obtain 



= 2r'^h'[ + 2(n + z)rh\ + zr {g[ - f[) + 2{n - l)zgi, 

= 2{z- l)rh[ + {n- l)rg[ + [z{z - 1) + n{n - 1)] gi - (z - l)(n + ^ - l)(/i - 2hi), 
0^2{z- l)rh[ + {n- l)rg[ + [z{z - 1) + n{n - 1) + 2{n - l){z - 1)B] gi 
-{z-l){z-n + l){fi-2h,) 



where B — {l^ — 4:a2p + 9as)/L^. Note that all the parameters of Lovelock gravity 
are in B, with B = 1 in Einstein gravity. The solution of the above equations may 
be written as 



7 = {(17 - 8B)z'^ - 2(3n + 9 - 8B)z + + 6n + 1 - 8B} ' , 
Fi = -2(^-1) (n-z -l))C-\ 
F2 = {J^i - H {8zK \{z -l)B + 2n + z- S]}'^ , 
F3 = {Ti + J2) {8zJC [{z -l)B + 2n + z- S]}"' , 
Gi = 2{z-l){n + z-l)IC'\ 
G2 = {Gi + Q2) {{z -l)B + 2n + z- 3]}^' , 
Gs = {Qi - G2) {8z/C \{z -l)B + 2n + z- 3]}"' , 

}C = {z-l){n + z-l)B + z{z - 1) + n(n - 1), 
J^i = 8{z - l)[{z - l)(n + z-l)B + z{z - 1) + n(n - 1)] 

X [{z - l)(n + 3^ - 3)B - 2z^ + (n + 3);^ + n{n - 2) - 1], 
jF2 = 7[n - 1 + - 1)B] |8(1 + B){z- if 

+{17n - 9 + 8B){z - if + 2{n + 8){n - l){z - 1) + n^(n - 1) - (n - 1)7^}, 

Qi = 8{z - l)[2{z - 1)B - 3^ + 3n - l][{z - l)(n + z - 1)B + z{z - 1) + n{n - 1)], 
02 = {8(1 + B){z- If + (17n ~9 + 8B){z- if 

+2{n + 8){n- l){z - 1) + n^{n - 1)|7 - (n - 1)7^ 



/(r) = l + £/i(r), 
g(r) = l + egi{r), 
h{r) = 1 + £/ii(r). 



Mr) 
9i{r) 



where Ci, C2 and C3 are integration constants and 
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Note that the functions fi{r), gi{r) and hi{r) reduce to those given in ref. [Q for 
«2 = «3 = (i3 = 1) and n = 3. For arbitrary values of a2 and da, 7 is not an 
integer, and therefore the only integer power of r is r"""^"*"^. However in Einstein 
gravity 7 can be an integer. This occurs for z = 2 with either n = 3 or n = 6. For 
z = 2 and n = 3, 7 = 4 and C2 = 0; the other two powers are the same and equal 
to —4. Of course, the next order leading term is 7'2(-"-2+i) = ^-8 -^yhich is a second 
order perturbation in e. For z = 2 and n = 6, 7 = 5. The largest power of r in the 
large r expansion is r~^. For other z or n, even in Einstein gravity, the only integer 
power of r at large r for = to first order in e (powers larger than 2{—n — z + 1)) 



IS 1^ ^ ^ ^ ^ 



We next obtain the coefficients of integer powers of r for k = ±1 up to the first 
order in e. The functions f{r), g{r) and h{r) at large r up to the first order in e may 
be written as 



2(n+2)-3 



/(r) = l+ Yl 



i=l ^ 
2(n+2)-3, 



9[r) 



1+ E - 

i=l 
2(ra+2)-3 



h{r) = l+ Yl 



i=l 



First, we consider the problem in Einstein gravity. For n = 3 and z = 2, all the 
coefficients of integer powers of r up to the first order [n > —7) are proportional to k 
except the coefficient of r~^. These coefficients are given in[|13]. For z = 2 and n = 6, 
the first non-vanishing term is and all the odd and even powers of r are present. 
It is a matter of calculation to show that the coefficients of first three powers (r~^, 
and r~^) are 

ai = 3ci, bi = y, 

02 = + ^kP, 62 = -^cj + ^kl^, C2 = -^cj + ^kf, 

1 3 136, ,2 J, 73 8 , ,2 29 3 44, ,2 

as = CT H kl Ci, 03 = -Ci H kl Ci, C3 = — C] H /c< Ci, 

^ 27 ^ 45 ' 9 ^ 15 ^' ^ 54 ^ 45 

where Ci is an arbitrary coefficient. For other values of n with z = 2 m Einstein 
gravity, the nonzero coefficients are those of r~^, r~"~^+-^, r~"~^~^,... 7--"-z+i-2m^ 
where m is the largest integer less than (n + z — l)/2 up to the first order terms in e 
{i < 2(n + 2; — 1)). Of course, the next power of r is 7--2't--22+2^ which is a second order 
term in our analysis and is nonzero even for k = case. We note that for 2; = n — 1 
in {n + l)-dimensional Einstein gravity, some logarithmic terms will appear. 



- 22 - 



For arbitrary values of the Lovelock coefficients, all the odd powers of 1/r vanish 
at large r and the coefficients of are 

a2 = k—{3a3[2z^ - 2{n - 5)z'^ + 3(n - 4)(n - 5)z - (2n^ - lln + 16)] 
B 

-2fa2[z^ + (n - 5)^2 + {n^ -7n+ ISfz - {2n^ - lOn + 13)] + (n - 2)[(n - 3)z - 2n + 5]/^}, 

62 = A;— {3q;3[2z^ + 2(n - 5)^^ + (n^ - lOn + 22)^ - 2(n - 3)^] + (n - 2)[(n - 4)z - 2(n - 3)]Z^ 

-2Z2q;2[^^ + (n - 5)^2 + (n - 3)(n - 5))^ - (2n^ - lln + 15)]}, 
C2 = A;— {3a3[2^3 + 2(n - 4)^^ ^ (^^2 _ 9^ + ^g)^ _ 2(^ - 2)(n - 3)] 

2x5 

-1l^dL2\z^ - (n - 4)^2 + {n^ -7n + n)z - (2n^ -9n + 10)] - (n - 2)[(n - 3)z - 2{n - 2)]Z^}, 
where 

B = (z+n-3) {Sasiz -2){z + n-3)- 2d2p[z{n - 3) - 2n + 5] - [z'^ - {n - l)z + 2{n - 2)] . 

The higher order coefficients with even power of r (r~^, r~^,..) are present and can be 
easily calculated but we shall not write them here as they are quite lengthy. In this 
case all a„'s are proportional to k. Thus, this feature is again a difference between 
Einstein and Lovelock gravity. That is, in Lovelock gravity all the even powers of 
r up to r"""^"*"^ are present in the expansion of the functions at large r, while in 
Einstein gravity this does not occur for 2; = 2. 
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